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^ . Abstract 

. The solution of QCD equations for generating functions of mul- 

ly-^ I tiplicity distributions reveals new peculiar features of cumulant mo- 

0^ ' ments oscillating as functions of their rank. This prediction is sup- 

r~| ' ported by experimental data on e^e~ ,hh, AA collisions. Evolution of 

P^, the moments at smaller phase space bins leads to intermittency and 

Q^l fractality. The experimentally defined truncated generating functions 

^ ■ possess zeros in the complex plane of an auxiliary variable recalling 

Lee- Yang zeros in statistical mechanics. 
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b : 1 Introduction 



For a long time, the phenomenological approach dominated in description 
of multiplicity distributions in multiparticle production. The very first at- 
tempts to apply QCD formalism to the problem failed because in the simplest 
double-logarithmic approximation it predicts an extremely wide shape of the 
distribution that contradicts to experimental data. Only recently it became 
possible to get exact solutions of QCD equations which revealed much nar- 
rower shapes and such a novel feature of cumulant moments as their oscil- 
lations at higher ranks. These moments are extremely sensitive to the tiny 
details of the distribution. Surprisingly enough, those QCD predictions for 
parton distributions have been supported by experimental data for hadrons. 



QCD is also successful in qualitative description of evolution of these 
distributions with decreasing phase space bins which gives rise to notions of 
intermittency and fractality. However, there are some new problems with 
locations of the minimum of cumulants at small bins. 

The experimentally defined truncated generating functions possess an 
intriguing pattern of zeros in the complex plane of an auxiliary variable. It 
recalls the pattern of Lee- Yang zeros of the grand canonical partition function 
in the complex fugacity plane related to phase transition. 

Before demonstrating all these peculiarities let us define the multiplicity 
distribution 

oo 

Pn = CFnl CFn-, (l) 
n=0 

where is the cross section of n-particle production processes, and the 
generating function 

oo 

G{z) = Y. Pr.il + zr. (2) 

ra=0 

The (normalized) factorial and cumulant moments of the P„ distribution are 
^ _ Y.nPnn{n-l)...{n-q + l) _ 1 d^G'(^) , 

where (n) = Y^nPnTt is the average multiphcity. They describe full and 
genuine g-particlc correlations, correspondingly. Let us point out here that 
the moments arc defined by the derivatives at the origin and are very sensitive 
to any nearby singularity of the generating function. 

In practice, one deals with distribution truncated due to finiteness of the 
available phase space and the summation in all formulae above is cut off 
at some finite value of n = N„iax which depends on the phase space region 
chosen, and increases with its increase. It is a polynomial of the power N^^x 
and has Nmax zeros in the complex 2;-plane. 

To shorten the presentation, I omit here all the details of calculations and 
references to original papers. The reader can find them in my review paper 
in Physics-Uspckhi 37 (1994) 715. Main qualitative results are described 
and demonstrated in Figures in the subsequent three Sections. Their physics 
implications are discussed in the last Section. 



2 Oscillations of cumulants of multiplicity dis- 
tributions in QCD 

First, let us consider QCD without quarks, i.e. gluodynamics. The gener- 
ating function of the gluon multiphcity distribution in the full phase-space 
volume satisfies the equation 

9G{z Y) ^ clxK{x)-f^[G{z,Y + lnx)G{z,Y + ln{l-x))-G{z,Y)]. (5) 
oY Jo 

Here Y = \n{p9/Qo),p is the initial momentum, 6 is the angular width of the 
gluon jet considered, p9 = Q where Q is the jet virtuality, Qq =const. 



7? = (6) 



TT 



^5 is the running coupling constant, and the kernel of the equation is 

K{x) = ^- {l-x)[2-x{l-x)]. (7) 

It is the non-linear integro-differential equation with shifted arguments in 
the non-linear part which take into account the conservation laws, and with 
the initial condition 

G{z,Y = 0) = l + z, (8) 

and the normalization 

G{z = 0,Y) = l. (9) 

The condition (j^) normalizes the total probability to 1, and the condition 
(H) declares that there is a single particle at the very initial stage. 

After Taylor series expansion at large enough Y and differentiation in eq. 
(H), one gets the differential equation 

(InG(F))" = 7o'[G'(r) - 1 - 2hG'{Y) + h^G'^Y)], (10) 

where hi = 11/24; /i2 = (67 — 67r^)/36 ~ 0.216, and higher order terms have 
been omitted. 

Leaving two terms on the right-hand side, one gets the well-known equa- 
tion of the double-logarithmic approximation which takes into account the 



most singular components. The next term, with hi, corresponds to the mod- 
ified leading-logarithm approximation, and the term with /12 deals with next- 
to-leading corrections. 

The straightforward solution of this equation looks very problematic. 
However, it is very simple for the moments of the distribution because G{z) 
and lnG{z) are the generating functions of Fg and Kg, correspondingly, ac- 
cording to (^, @. Using this fact, one gets the solution which looks like 

^ Kg ^2[i_2/iig7 + /^2(gV + gy)] .... 
Fg- + qy ' 

where the anomalous dimension 7 is related to 70 by 

7 = 70 - lha'o + ^(4/^2 - hlhl + 0(7o'). (12) 

The formula ([llD shows how the ratio Hg behaves in different approxi- 
mations. In double-log approximation when hi = h2 = 0, it monotonously 
decreases as that corresponds to the negative binomial law with its pa- 
rameter k = 2 i.e. to very wide distribution. In modified-log approximation 
(/;,2 = 0) it acquires a negative minimum at 

qmin = r— + l + 0(7o) ~ 5 (13) 

and approaches asymptotically at large ranks q the abscissa axis from below. 
In the next approximation given by (|Tl|) it preserves the minimum location 
but approaches a positive constant crossing the abscissa axis. In ever higher 
orders it reveals the quasi-oscillatory behavior about this axis. This predic- 
tion of the minimum at g ^ 5 and subsequent specific oscillations is the main 
theoretical outcome. 

It is interesting to note that the equation (^ can be solved exactly in 
the case of fixed coupling constant. All the above qualitative features are 
noticeable here as well. 

While the above results are valid for gluon distributions in gluon jets 
(and pertain to QCD with quarks taken into account), the similar qualitative 
features characterize the multiplicity distributions of hadrons in high energy 
reactions initiated by various particles. As an example, I show in Fig.l the 
ratio Hg as a function of q in the e"'"e~ data of DELPHI collaboration at 
91 GeV, where the oscillations and the location of minima are of a special 
interest. 



3 Evolution of distributions with decreasing 
phase-space volume — intermittency and frac- 
tality 

The multiphcity distributions can be measured not only in the total phase 
space (as has been discussed above for very large phase-space volumes) but in 
any part of it. For the homogeneous distribution of particles within the vol- 
ume, the average multiplicity is proportional to the volume and decreases for 
small volumes but the fluctuations increase. The most interesting problem 
here is the law governing the growth of fluctuations and its possible depar- 
ture from a purely statistical behavior related to the decrease of the average 
multiplicity. Such a variation has to be connected with the dynamics of the 
interactions. In particular, it has been proposed to look for the power-law 
behavior of the factorial moments for small rapidity intervals 5y 

(X (0(g) >0) (51/^0), (14) 

inspired by the idea of intermittency in turbulence. In the case of statistical 
fluctuations with purely Poissonian behavior, the intermittency indices 0(g) 
are identically equal to zero. 

Experimental data on various processes in a wide energy range support 
this idea (e.g., see Fig.2), and QCD provides a good basis for its explanation 
as a result of parton showers. The generating function technique is not 
applicable here, and one should consider the Feynman graphs of evolution 
of a jet with its subjet hitting the phase-space window under consideration 
(see the abovementioned review). 

At moderately small rapidity windows, one can get in the double-log 
approximation the power-law behavior with 

0(g) = D(g-l)-^^7o. (15) 

The running property of QCD coupling constant is not important in that re- 
gion. This property becomes noticeable at ever smaller windows when (e.g., 
at g=2) \n.5yQ/5y > a^^, and leads to smaller numerical values of 0(g) com- 
pared to ([T5|). The general trends in this region decline somewhat from the 
simple power law (|I^) due to logarithmic corrections. Qualitatively, these 



predictions correspond to experimental findings at relatively small ranks q 
where the steep increase in the region oi 5y > 1 on the log-log plot of the 
dependence (O) is replaced by slower one at smaller intervals 5y (see Fig. 2). 



The transition point between the two regimes depends on the rank in quali- 
tative agreement with QCD predictions also. Namely, the transition happens 
at smaller bins for higher ranks. These findings can be interpreted as an indi- 
cation on fractal structure of particle distributions within the available phase 
space. When interpreted in terms of fluctuations, they show that the fluctu- 
ations become stronger in small phase-space regions in a definite power-like 
manner and, surely, exceed trivial statistical fluctuations. 

Let us turn now to the g-behavior of moments at small bins. The phe- 
nomenon of the oscillations of cumulants discussed above reveals itself here 
as well if one goes beyond the double-log approximation of (|15]). In terms of 
factorial moments, it means the non-monotonous behavior of the intermit- 
tency indices as functions of q. (Compare it to the steady increase with q at 
g > 1 given by (|15|).) It gives rise to the negative values of Kq and Hq. 

The fate of the first minimum can be easily guessed from the formula 
(p!3D. For large enough virtualities (i.e. small 70), the minimum location 
moves to higher values of rank q for jets with larger virtuality Q since the 
QCD coupling constant is running as ln~^ Q. Therefore, the predicted shift 
of the minimum is 

In^/^Q. (16) 

It follows that qmin nioves to higher ranks at higher energies because more 
massive jets become available. Another corollary is that it should shift to 
smaller values of q for smaller bins at fixed energy. 

While former statement finds some support in experiment, the second one 
does not look to be true as shown in Fig.3. On the contrary, the minimum 
appears at higher ranks for smaller bins. There is no solution of this problem 
yet but it should be ascribed to the higher-order effects. Actually, one can 
guess that the higher order terms shown as 0(70) in (|l^) become so important 
at small bins that they overpower the weak Q-dependence of 75"^ in the first 
term of (pIS]). It is important to stress here that at large rapidity intervals the 
modified leading-log term with h2 does not influence the value of qmin, and 
only increases the value of Hq by /i27o- Thus, the next- to- leading corrections 
should be in charge of the additional shift of qmin, and, therefore, small bins 
help us look into higher orders of QCD. 



4 Zeros of truncated generating functions 



There is another fascinating feature of multiphcity distributions - it happens 
that zeros of the truncated generating function form a spectacular pattern 
in the complex plane of the variable z. Namely, they seem to lie close to a 
single circle. At enlarged values of Nmax they move closer to the real axis 
pinching it at some positive value of z. It is demonstrated in Fig. 4 for UA5 
data on pp interactions at 200 and 900 GeV for various rapidity windows. 

No QCD interpretation of the fact exists because it is hard to exploit 
the finite cut-off in analytic calculations. The interest to it stems from the 
analogy to the locations of zeros of the grand canonical partition function as 
described by Lee and Yang who related them to possible phase transitions 
in statistical mechanics. In that case, z variable plays the role of fugacity, 
and pinching of the real axis implies existence of two phases in the system 
considered. 

In particle physics, it shows up the location of the singularity of the 
generating function i.e. the number of zeros of truncated generating functions 
increases and they tend to move to the singularity point when Nraax ~^ 
oo. Since it happens to lie close to the origin, it drastically influences the 
behavior of moments (see (^), (p), and, therefore, determines the shape of 
the distribution. The study of the singularities is at the very early stage now, 
and one can only say that the singularity is positioned closer to the origin 
in nucleus-nucleus collisions and it is farthest in e~^e~ that appeals to our 
intuitive guess. 

5 Discussion and conclusions 

Let us discuss the implications of the above findings. The QCD prediction 
of quasioscillating behavior of cumulant moments of the multiplicity distri- 
butions reveals the tiny features which were overlooked by simpleminded fits 
of the negative binomial distribution. The truncated negative binomial dis- 
tribution has the oscillating cumulants as well but the oscillations should die 
out asymptotically while they persist in QCD. 

It demonstrates that the QCD distribution belongs to the class of non- 
infinitely-divisible ones and shows that the Poissonian cluster models are 
ruled out by QCD since they lead to positive values of cumulants and can 



not reproduce their oscillations at asymptotically high energies. 

Moreover, this prediction is closely related to the existence of a new ex- 
pansion parameter in describing multiparticle production in QCD. In terms 
of moments, this parameter is equal to the product of the moment's rank 
to the anomalous dimension of QCD 7^. Let us recall that the similar pa- 
rameter appears when calculating the Feynman tree graphs and looks like 
a product of the number of final particles to QCD coupling constant nas- 
Both parameters become large if large number of particles is involved. 

The invalidity of simple perturbative approach in multiple production 
processes would ask for more convenient basis (compared to particle number 
representation) to be used. As one of possible examples of effectiveness of 
such an approach, we would mention lasers where coherent state basis is 
more suitable. That is why studies of coherent and squeezed states as well as 
statistical approach and search for collective effects, in general, are welcome. 
The recent findings of zeros of the truncated generating functions point out 
in that direction, and provoke speculations to their possible relation to the 
problem of phase transition in hadronic matter. 

The success of QCD in predicting the qualitative features of moments 
of multiplicity distributions both for large phase-space regions (oscillation of 
cumulants) and for small bins (intermittency indices, fractality) seems even 
more surprising if one recognizes that they are derived for parton distribu- 
tions while hadrons are observed in experiment. Nevertheless, these quali- 
tative features reveal themselves in hadron distributions as well. Why such 
a local parton-hadron duality persists even at higher moments is an open 
question. Another problem appears in such a tiny property as the evolution 
of the minimum location at smaller bins. Its solution can show a way to take 
into account higher order effects of QCD more properly. 

To conclude, we would hke to stress that recent theoretical studies of 
QCD multiplicity distributions predict rather exotic features of the moments 
of the distributions which find some support in experiment and, at the same 
time, provoke new problems to be solved. 
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FIGURE CAPTIONS 

Fig.l The ratio of cumulant to factorial moments Hq as the function of 
the rank q. The DELPHI-data on e"'"e~ at 91 GeV in the total phase space 
are shown by dots. The dashed line shows the fit by the negative binomial 
distribution with the parameters given by DELPHI. The sohd line is drawn 
just to guide the eye. 

Fig. 2 Factorial moments of order ^=2,. . .,4 for the all-charged, positives- 
only and negatives-only samples of NA22-data on Tz^p interactions at 250 
GcV/c (lab. system). 

Fig. 3 Behavior of the ratio of cumulant to factorial moments Hq as a 
function of the rank q for some single hemisphere multiplicity distributions 
from the DELPHI collaboration: (a) full hemisphere; (b) rapidity window 
\y\ < 1-5; (c) rapidity window \y\ < 1.0. The minimum shifts to higher q for 
smaller bins. 

Fig. 4 The locations of zeros of the truncated generating function for UA5- 
data on pp interactions at 200 and 900 GeV (c.m.s.) in various rapidity 
windows. The upper halfplane oi z — x -\- iy is shown only because of the 
up-down symmetry. (For further experimental information see the talk of G. 
Gianini at Multiparticle Dynamics-95). 
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